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I.  INTRODUCTION 


Not  only  is  the  physical  interpretation  of  the  radiation  of  the  acoustic  waves  in 
more  than  one  dimension  complex,  but  the  effects  of  varying  frequency  and  shading 
are  hard  to  visualize  when  considering  directionality  and  frequency  properties.  As 
Wood  implied  in  his  book,  A  Textbook  of  Sound  [Ref.  1:  p.449j,  the  directional  factor 
is  very  important  in  many  acoustic  applications  such  as  speech,  hearing  and  array 
design. 

The  simplest  array  one  can  consider  is  the  beam  pattern  of  a  dipole  radiating  in  two 
dimensions.  References  to  the  case  of  the  simple  dipole  can  be  found  in  virtually  every 
textbook  in  acoustics.  For  a  fairly  recent  example,  see  Malecki's  book  [Ref.  2:  p.l99j. 
Currently,  the  mathematical  modeling  of  the  simple  acoustic  arrays  has  been  enhanced 
with  the  advent  of  the  sophisticated  computer  software.  Specifically,  Ziomek 
[Ref.  3:  p.  1 38]  and,  Embleton  and  Thiessen  [Ref.  4:  p.  1 124]  have  investigated  N  point 
source  arrays,  considering  both  amplitude  and  phase  weighting.  Although  the 
mathematics  and  computer  hardware  exist,  the  visualization  of  the  three-dimensional 
aspects  of  these  arrays  is  difficult.  By  use  of  the  computer  graphics  software  the  ability 
to  interpret  the  significance  of  shading  and  phasing  can  be  fully  developed  [Refs.  5,6]. 

To  date,  visualization  of  the  three-dimensional  beam  pattern  has  been  needed  to 
understand  the  directional  factor.  Numerical  computer  output  was  essential  to  draw 
three-dimensional  graphs.  By  changing  two  parameters  (in  three-dimensions,  those 
could  be  x  and  y  or  0  and  (p,  etc.),  different  generated  values  of  the  beam  pattern,  F, 
are  obtained. 

One  computer  program  to  visualize  the  beam  pattern  in  three-dimensions  (P3D)  and 
four  computer  programs  to  read  values  in  detail  (2DICIR,  2DPCIR,  2DIFIX, 
2DPFIX)  were  written.  These  computer  models  predict  the  directional  factor  and 
overall  beam  pattern  with  or  without  differences  in  the  source  strengths,  inter-source 
distances,  and  phases  of  the  elements  of  the  arrays. 

Typical  geometry  for  the  beam  pattern  prediction  is  shown  in  Figure  2.3  on  page 

16. 

The  following  definitions  apply  and  will  be  used  throughout: 

1  *  source  strength  of  the  first  source  (source  1) 
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A  »  source  strength  of  source  2 

B  -  source  strength  of  source  3 

C  •  source  strength  of  source  4 

D1  *  distance  from  the  geometrical  center  to  source  1 

D2  “  distance  from  the  geometrical  center  to  source  2 

D3  *  distance  from  the  geometrical  center  to  source  3 

D4  =  distance  from  the  geometrical  center  to  source  4 

<p2  =  phase  difference  between  source  2  and  source  1 

<p3  *  phase  difference  between  source  3  and  source  l 

<p4  *  phase  difference  between  source  4  and  source  1 

0  =  angle  measured  from  the  z  axis  toward  the  x-y  plane 

<I>  =  angle  measured  counterclockwise  in  the  x-y  plane  starting  from  the  x  axis 

k  -  wave  number 

w  =  angular  frequency 

R  *  distance  from  the  geometrical  center  to  a  certain  point  in  space 
H  =  directional  factor 

F  *  beam  pattern  =  201og10|H|  or  101og10H2  expressed  in  dB 

All  distances  are  measured  from  the  origin  of  the  coordinate  system  to  the  sources 
along  the  x  and  y  axes.  Phase  differences  are  based  upon  the  phase  of  the  first  source, 
which  is  assumed  to  be  zero. 

The  goal  of  this  thesis  is  to  develop  a  convenient  methodology  for  the 
investigation  of  the  functional  dependency  of  the  directional  factor  on  the  source 
amplitudes  (A,  B  and  C),  distances  from  the  origin  (Dj,  D2,  D3  and  D4),  phase 
differences  (<p2,  <p3  and  <p4),  and  0  and  <I>  as  demonstrated  by  the  five  programs:  P3D, 
2DICIR,  2DPCIR,  2D1FIX,  2DPFIX. 


II.  THEORY 


A  spatial  configuration  of  simple  sources,  either  discrete  or  distributed,  each  with 
its  own  complex  source  strength,  can  be  used  to  represent  the  quadruplet  source  case. 
The  pressure  at  a  field  point  is  the  sum  of  pressures  produced  by  the  individual  sources 
[Ref.  7:  p.  169].  Before  we  develop  the  quadruplet  source  case,  the  acoustic  doublet 
case  needs  to  be  described  to  provide  a  background  for  the  more  complicated 
quadruplet  case. 

A.  DOUBLET  SOURCES 

A  doublet  source  consists  of  two  point  sources  of  strengths  of  1  and  A,  separated 
from  the  origin,  by  distances  Dj  and  D2  and  vibrating  at  the  same  frequency  but  with 
a  phase  difference  <p. 


From  the  geometry  shown  in  Figure  2.1,  we  have 


cos0  *  8j/ Dj, 
sin0  -  Cj/Dp 


Tj  =  V(R  •  £j)2  +  V 


-*/(R  •  DjSinQ)2  +  (DjCOsO)*, 

*  V(1  -  DjSinQ/R)2  •+•  ^cosO/R)2. 


R 

sin0  *  e2/D2, 
cos0  =  6,/D,f 


r,  =  V(R  +  c?)2  +  6, 2 


V(R  •**  D.,sin0)2  +  (D,cos0)2, 


^-.va 


D,sin0.'R)2  +  (8,cos0,'R)2. 


The  pressure  at  point  (R,0)  due  to  source  1  is 

P.  =  —  ei(w  -  ^i) 
ri 

where  w  is  the  angular  frequency.  That  due  to  source  2  is 


(2.1) 

(2.2) 


(2.3) 

(2.4) 

(2.5) 

(2.6) 


(2.7) 

(2.S) 


(2.9) 


P2 


—  ei(wt  -  kr2  -  (p) 
r2 


(2.10) 


The  total  acoustic  pressure  at  point  (R,0)  is  then  the  sum  of  Pj  and  P2, 


p  =  _L  ei(;vt  -  kR){  —  eik(R  ‘  ri) 


e"ilk( r2  -  R)  +  <Pl  } 

r2 


(2.11) 


where  R  is  the  distance  from  the  field  point  to  the  origin.  From  (2.4)  and  (2.8),  we  get 


k(R  -  Tj)  =  kR(  1  -  rj/R) 


.a  f.«  r.« 


>t 


-  kR[l  -  V(1  -  DjSine/R)2  +  (DjCosQ/R)2  ],  (2.12) 

k(r2  -  R)  -  kR(r2/R  -  1) 

-  kRt^DjSine/R  +  1)^  +  (D2cos0/R)i  -  1  ].  (2.13) 

Therefore,  equation  (2.11)  can  be  rewritten  as 


p  =  _L  ej(wt  -  kR)j  Ji  ejkR[  1  -  V(i  -  DjSinS/R)2  +  (DjCosG/R)2  ] 

R  rl 


_^e-j[kR[V(D2sin0/R  +  l)2  +  (D2cos0/R)2  •  1  )  +  <p]  j  ^2 


14) 


In  the  most  frequently  encountered  cases,  the  observation  of  the  pressure  is  made 
at  distances  greater  than  the  separation  of  the  sources.  Therefore,  it  will  be  useful  to 
derive  the  far  field  (R>  >  (Dj  +  D2))  pressure  equation  from  the  above  equation  (2.14). 
From  the  geometry  shown  in  Figure  2.2,  we  have 


sin0  =  Ej/Dj  =  £2/D,. 


(2.15) 


Therefore,  (R  -  rt)  and  (r2  -  R)  can  be  replaced  with  DjSin©  and  D2sin0.  Also  we  have 


(2.16) 


so  that  the  pressure  equation  becomes 


p  „  _L  *  kR^gjkDjSinO  +  Ae-j(kD2sin0  +  <p)j 


(2.17) 


IPI  =  P(R)H(0) 


(2.18) 


where  P(R)  =  (1+A)/R. 
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B.  QUADRUPLET  SOURCES 

Similarly,  a  quadruplet  source  case  can  be  developed  with  four  sources  of 
strengths  1,  A,  B  and  C,  separated  from  the  origin,  by  distances  DJt  D2,  D3  and  D4 


R 


t: 


R 


Figure  2.3  Geometry  of  quadruplet. 
From  the  geometry  shown  in  Figure  2.3,  we  have 


x  *  Rsin8cos<I>, 
y  *»  RsinGsinO, 
z  «  Rcos8, 


R  -  Jtc  +  y2  +  z2 
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(2.21) 

(2.22) 

(2.23) 


(2.24) 


where  x,  y,  z  are  rectangular  coordinates  and  R  is  the  distance  from  the  field  point  to 
the  origin. 

And  then 


rt  ■  V(y  +  Dj)2  +  x2  +  z2. 

(2.25) 

By  Taylor's  approximation 

h  -  R  +  Djy/R  +  D^R. 

(2.26) 

rj  -  R  »  Djy/R  -  DjSinOsinO, 

(2.27) 

r2  *  V(y  -  D2)2  +  x2  ■+•  z2  -  R  -  D2y/R, 

(2.28) 

R  -  2  "  D2y/R  ■  D2sin0sin<l> 

(2.29) 

where  R  is  much  greater  than  D1(  and  D2. 

Using  the  same  approach 

(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 

where  R  is  much  greater  than  D3,  and  D4. 

Therefore,  the  pressure  at  (R,  0,  <p)  is 

p  =  _L  ej(wt  -  kR  -  kDjSin6sin<J>)  +.  A  ej(wt  -  kR  +  kD2sin0sin<t>  -  <p,) 

R  R 

+  £  ej(wt  -  kR  +  kD3sin0cos<l>  -  <p3) 

R 

+  £  ej(wt  -  kR  -  kD4sin0cos<t>  -  <p4) 

R 


r3  *  V(x  *  D3)2  +  y2  +  z2  -  VR2  -  2D3x  +  D32 
-  RVl  -<2D3x/R2-  D^.'R2). 

r3  -  R(1  -  (1/2X2D3x/R2  -  D32/R2)]  -  R  -  D3x/R  +  D32'2R, 
R  -  r3  *  D3x/R  ■  D3sin0cosO 

r4  -  V(x  +  D4)2  ~  y2  +  z2  -  R  +  D4  x;R, 
r4  -  R  *  D4sin0cos<P 


(2.35) 
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The  directivity  factor, 


H(0  O)  -  _ - _ I  e'i(kD l  sm©sin<I>>  +  Aej(kD2sin0sinO  -  <p2) 

(1  +  A  +  B  +  C) 

+  Bej(kD3sin0cosO  *  <P3)  +  Ce'i^kD4sin0cos<I>  +  9*)  |. 
Expanding  the  complex  terms  yields 


(2.36) 


H(0,O) 


(IVa  +  bVc)  ’^cos^k^ls^n0s^n<^^  +  Acos(kD2sin0sinO  -  <p2) 

+  Bcos(kD3sin0cosO  -  <p3)  +  Ccos(kD4sin0cos<I>  +  q>4)]2 
+  (su^kDjSinBsinO)  -  Asin(kD2sin0sinO  -  <p2) 

-  Bsin(kD3sin0cosO  -  <p3)  +  Csin(kD4sin0cosO  +  <p4)]2} </ 2 . 


(2.37) 


This  is  the  expression  used  for  the  computer  program. 

C.  GRAPHICAL  PRESENTATION  OF  THE  DIRECTIVITY 

1.  Two-Dimensional  Rectangular  and  Polar  Plots 

From  equation  (2.37),  H  (0,<P)  versus  0,  and  H  (B,^)  versus  4>  can  be  plotted 
if  either  <D  or  0  respectively  is  given  constant  value.  In  a  rectangular  plot,  the 
directivity,  F,  which  is  equal  to  20log10|H  (0,<I>)  |  is  plotted  on  one  axis  and  0  or  <P  on 
the  other.  An  example  of  this  is  seen  in  Figure  3.3.  In  a  polar  plot,  the  directional 
factor  or  beam  pattern  is  a  magnitude  of  the  radial  distance  from  the  center  of  the  plot 
while  the  appropriate  angle  is  measured  clockwise  around  the  origin.  An  example  of 
this  is  seen  in  Figure  3.4. 

2.  Three-Dimensional  Cartesian  Plot 

From  equation  (2.37),  H  can  be  calculated  and  plotted  as  a  function  of  both  0 
and  <D.  This  provides  the  three-dimensional  surface  with  its  elevation  equal  to 
201og10|H  (0,0)1  for  each  0-0  pair.  An  example  of  this  is  Figure  3.7.  This  form  of 
presentation  is  difficult  to  interpret  but  can  be  a  very  useful  analytical  tool.  Further 
explanations  are  given  in  the  next  chapter. 

3.  Contour  Plot 

A  contour  plot  gives  lines  of  constant  H  (0,0).  As  an  example,  let  kDj  * 
kD2  =  kD3  =  kD4  *  kD,  <p2  *  0,  <p3  =*  tp4  »  n  and  A  =  B  =  C  =  1.  Equation 
(2.37)  then  can  be  simplified  to 


III.  PROGRAM  DEVELOPMENT 


In  analyzing  the  directional  and  frequency  properties  of  shaded  and  phased 
simple  arrays  several  types  of  display  methods  were  needed.  At  the  first  stage  of 
developing  the  computer  program,  a  HEWLETT-PACKARD  86B  computer  and  a 
HP709QA  MEASUREMENT  PLOTTING  SYSTEM  were  used.  By  using  this 
equipment  we  were  able  to  get  two-dimensional  plots  as  shown  in  Figures  3.1  and  3.2. 
However,  it  took  much  time  in  executing  the  program  and  also  we  could  not  do  three- 
dimensional  plottings. 


Figure  3.1  2-D  rectangular  plot  (HP  86B). 
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Figure  3.2  2-D  circular  plot  (HP  86B). 

To  get  a  complete  understanding  of  the  beam  pattern,  three-dimensional  displays 
viewed  from  various  angles  were  essential.  Consequently,  alternative  methods  of 
display  were  programmed  in  FORTRAN  and  executed  on  the  IBM  3033.  There  are 
several  reasons  why  this  proved  to  be  the  most  desirable  method.  Primarily,  a  three- 
dimensional  graphics  software  package  (DISSPLA)  is  available  for  use  on  the  IBM 
3033.  Secondly,  FORTRAN  generated  the  numencal  output  with  greatly  increased 
efficiency  when  using  the  EXTENDED  FORTRAN  COMPILER  l Ref.  8:  p.lS], 
Finally,  the  IBM  3033  has  many  more  output  devices  such  as  VERSTEC.  TEK61S  and 
SHERPA  to  name  a  few.  Therefore,  the  programs  used  in  the  HEWLETT- 
PACKARD  86B  were  rewritten  in  FORTRAN.  The  translation  from  BASIC  to 
FORTRAN  required  some  changes  in  variable  names  along  with  the  use  of  some 
convenient  function  calls.  For  reference,  some  of  the  original  BASIC  programs  are 
attached  in  Appendix  A. 
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Figure  3.3  2-D  rectangular  plot  (DISSPLA). 

In  the  trial  of  analyzing  only  with  numerical  outputs,  BASIC  had  become 
ineffective  in  recognizing  details  due  to  the  very  large  number  of  computed  values.  At 
this  stage  the  graphic  package  DISSPLA,  available  on  the  IBM  3033  computer,  was 
used  in  conjunction  with  these  programs.  The  DISSPLA  package  is  a  library  of 
FORTRAN  subroutines  that  facilitate  data  plotting.  It  is  also  device  independent, 
which  means  that  it  does  not  rely  upon  features  particular  to  any  type  of  graphical 
device  and  has  its  own  extensive  symbol  and  character  generation  routines.  By  being 
device  independent,  programs  of  the  HEWLETT-PACKARD  86B  could  be  translated 
without  difficulty.  Therefore,  DISSPLA  was  used  to  produce  both  two-dimensional 
and  three-dimensional  graphs.  The  data  used  to  compute  graphs  could  be  generated 
within  the  FORTRAN  program  or  read  from  a  data  file.  The  former  method  was  used 
here  because  data  files  took  a  great  amount  of  file  space  due  to  the  large  volume  of 
values. 
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THEEA=35  DEGREE 
K=105 

Dl=05  02=0.5 
D3-0.5  D4-05 
P2-0  P3-0  P4»  0 
A-L  B»L  OL 


Figure  3.4  2-D  circular  plot  (D1SSPLA). 

Two-dimensional  graphs  were  initially  tried.  By  holding  one  parameter  constant 
and  with  the  use  of  some  array  variables,  plots  of  directivity  versus  angle  were 
obtained.  Examples  of  two  kinds  of  displays  are  shown  in  Figures  3.3  and  3.4. 
Programs  are  attached  in  Appendix  B. 

To  predict  the  real  shape  of  a  directivity  pattern  in  three-dimensions,  three- 
dimensional  graphics  had  to  be  developed.  In  trying  to  get  some  satisfying  three- 
dimensional  plotting,  all  methods  in  DISSPLA  were  tried.  They  are  SURFUN,  and 
SURMAT  for  regular  matrix,  scattered  points,  and  three-dimensional  vector  drawing. 

SURFUN  is  a  subroutine  that  draws  the  surface  from  points  in  the  z  direction  in 
conjunction  with  a  grid  in  the  x-y  plane  on  the  condition  that  z“f{x,y)  is  given.  One 
constraint  that  was  faced  in  applying  this  function  was  that  there  was  more  than  one  z 
value  for  a  specific  x-y  coordinate. 
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Figure  3.5  3-D  plot(SURMAT  for  scattered  points). 


The  next  trial  was  SURMAT  for  a  regular  matrix  (for  example,  a  rectangular 
three-dimensional  lattice  of  points).  For  the  same  reason  as  in  applying  the  SURFUN 
subroutine,  this  also  turned  out  to  be  non-applicable  even  though  all  data  points  were 
translated  into  a  regular  matrix. 

The  next  option  tried  was  SURMAT  for  scattered  points  (for  example,  a  three- 
dimensional  irregular  cluster  of  points).  As  seen  in  Figure  3.5  it  also  was  unsuccessful 
due  to  more  than  one  point  on  one  specific  coordinate  of  the  x-y  plane.  The  program 
is  listed  in  Appendix  C  for  reference. 

The  last  option,  the  three-dimensional  vector  drawing  method  gave  us  a  more 
understandable  display  as  shown  in  Figure  3.6.  However,  it  did  not  show  the  beam 
pattern  clearly  because  all  the  lines  were  visible,  even  from  different  view  angles.  At 
this  stage  the  problem  was  the  removal  of  hidden  lines.  Unfortunately,  that  was  not 
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possible  because  the  starting  and  ending  points  of  the  hidden  lines  could  not  found. 
The  program  is  listed  in  Appendix  D. 
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K=10.55 
01=0.5  02=0.5 
D 3-0.5  D4-05 
VIEW  ANGLE”  45,  20 
P2-0  P3=180P4=180 
A=L  B=l.  C=L 


Figure  3.6  3-D  plot(vector  drawing  method). 


Consequently,  with  the  consideration  of  all  the  constraints  mentioned  above, 
another  type  of  format  was  made  to  see  the  details  of  the  beam  pattern.  Now,  a  plot 
was  made  using  the  function  SURMAT  about  the  0-<J>  plane  instead  of  the  x-y  plane. 
This  eliminated  the  previous  multiple  value  problem.  The  results  showed  the  most 
desirable  presentation  among  all  those  tried.  As  shown  in  Figure  3.7  it  gives  a  clearer 
appreciation  of  the  directivity  pattern  in  conjunction  with  0  and  <l>.  The  program  is 
listed  in  Appendix  E. 

We  can  relate  this  three-dimensional  graph  (Figure  3.7)  to  the  tw’o-dimensional 
graphs  (Figures  3.3  and  3.4)  for  interpretation.  As  shown  in  Figure  3.3  the  x-axis  is 
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the  angle  and  the  y-axis  is  the  directivity  which  is  expressed  in  dB.  Figure  3.4  is  the 
circular  graph  with  angle  working  counterclockwise  in  10°  increments  and  each  circle  is 
an  increment  of  10  dB  with  the  center  of  -  40  dB.  Figure  3.7  is  the  three-dimensional 
graph  from  the  view  point  of  20°  from  the  x-z  plane,  and  60°  from  the  x-y  plane.  Here, 
the  x-axis  is  <l>,  the  y-axis  is  0  and  the  z-axis  is  the  beam  pattern  (dB).  From  Figure 
3.3  four  major  points  can  be  found  at  approximately  50°,  130°,  230°  and  310°  of  0 
with  <J>  fixed  at  126°.  From  Figure  3.7  the  same  feature  can  be  determined  at  50°  of  0, 
126°  of  O  as  seen  in  Figure  3.3.  Since  Figure  3.3  was  made  by  holding  <t>  constant  at 
126°,  one  can  look  along  the  O  «  126°  line  in  Figure  3.7  and  note  the  null  at  0  = 
90°,  the  maximum  at  0  =  50°,  and  the  relative  maximum  at  0  =  0°.  Note  the  heavy 
lines  in  Figure  3.7. 
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K=105 

Dl=  05  D2=  0.5 
D3-05  D4-0.5 
VIEW  ANGLE— +20, +60 
P2-0  P3-0  P4-0 
A=l_  B=l.  C=l. 


Figure  3.7  3-D  plot(SURMAT  in  0-C>  plane). 
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From  the  circular  graph  (Figure  3.4)  four  major  points  can  be  found  approximatly  at 
40°,  140°,  220°  and  320°  of  <D  with  0  fixed  at  35°.  In  the  three-dimensional  graph  one 
of  these  major  points  is  seen  approximately  at  35°  of  0,  40°  of  0  and  the  other  major 
point  is  seen  approximately  at  35°  of  0,  140°  of  O.  Comparison  is  made  easier  when 
we  see  Figure  3.8  and  Figure  3.9  which  were  drawn  with  the  surfaces  cut  at  O  =  126° 
and  0  **  35°  respectively.  Note  the  shift  in  axes  between  the  two  figures.  All  these 
points  of  two-dimensional  graphs  correspond  to  those  of  the  three-dimensional  graph. 


Figure  3.8  3-D  plot(SL'RMAT  in  0-<I>  plane) 


Inversely,  every  point  of  the  three-dimensional  graph  can  be  read  in  two-dimensional 
graphs  to  see  detailed  values  whenever  needed. 

Finally,  a  trial  was  made  to  draw  a  constant  value  plot  of  the  beam  pattern.  The 
first  approach  was  to  derive  the  functions  of  <P  =  f  (0)  or  0  =  f  (<D)  and  plot  these  on 
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K=10.5 

Dl=05  D2=  0.5 
D3=  05  D4=  05 
VIEW  ANGLE=+20,+60 
P2=  0  P3=  0  P4=  0 
A=1  B=l.  0=1. 


Figure  3.9  3-D  plot(SURMAT  in  0-0  plane). 

the  0-<I>  plane.  However,  failing  to  derive  these  relationships  for  the  general  four 
element  case  made  this  method  of  making  a  contour  plot  of  a  fixed  F  value  impossible. 
Instead,  a  modified  version  of  the  existing  0*0  program  was  used.  It  sorted  on  a  small 
range  of  values  about  the  desired  value  of  F  to  find  the  corresponding  0-0  coordinates 
which  were  then  plotted  as  in  Figure  3.10.  This  type  of  plot  is  included  throughout  the 
study. 

With  these  forms  of  display  and  the  speed  of  the  FORTRAN  version  of  this 
program  we  can  now  use  this  method  to  study  the  directional  and  frequency  properties 
of  shaded  and  phased  simple  arrays. 
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IV.  COMPARISON  WITH  SAMPLE  CASES 


In  making  comparisons,  some  simple  parameters  were  used  to  make  the 
comparisons  easier.  There  are  three  cases  included  in  this  chapter.  Each  case  has  a 
unique  set  of  input  data,  and  computer  generated  and  experimental  plots.  In  all  cases, 
a  two  points  were  picked  and  hand  calculations  made  to  verify  the  computer  output. 
An  example  of  the  hand  calculations  is  given  for  the  first  case  only.  Comparisons  were 
also  made  between  hand  or  computer  generated  values  and  the  two-dimensional  polar 
plot,  between  the  theoretical  two-dimensional  plot  and  the  experimental  data  plot, 
between  the  two-dimensional  plot  and  the  three-dimensional  plot,  and  finally  between  a 
fixed  value  of  the  beam  pattern  (a  "track")  and  its  three-dimensional  plot. 

While  in  general,  0(0)  or  0(0)  for  a  constant  F  could  not  be  obtained 
analytically,  in  a  few  cases,  included  below,  it  was  possible  to  obtain  the  functional 
dependencies.  These  are  discussed  in  the  following  cases. 

A.  CASE  I 

When  kDj  =  kD2  =  kD3  =  kD4  =  5.275, 

<P2  =  0,  <p3  =  q>4  -  Tt. 

A  =  B  =  C  =  1 

equation  (2.37)  can  be  simplified  as 

H(0,O)  =  — —  |cos(5.275  sin0sinO)  -  cos(5.275  sin0cosO)|.  (4.1) 

Using  the  above  input  data  with  values  of  0  =  40°  and  O  =  0°  we  get  H(0,O)  = 
0.984566  and  F  (0,0)  =  -0.135104  from  the  computer  output  as  seen  in  Appendix  G. 
If  we  use  the  above  simplified  equation  and  hand-computation,  we  get  H(0,O)  = 
0.98457  and  F(0,O)  =  -0.135068  . 

In  the  same  way,  if  0  =  90°  and  0  =  0°  then  we  get  H  =  0.233302,  and  F  = 
-12.641622  from  the  computer  output  and  H  =  0.23330,  and  F  =  -12.6417  from  the 
hand-calculation.  These  hand-calculated  values  are  almost  the  same  the  the  computer 
output.  Comparisons  can  also  be  made  between  these  selected  values  and  the  two- 
dimensional  plot.  Figure  4.1a.  Additionally,  the  above  two-dimensional  plot 
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corresponds  well  in  shape  with  Figure  4.1b  which  is  from  LCDR  John  Butler's 
experiment  [Ref.  10]. 

If  we  compare  the  shape  of  the  two-dimensional  plot  (Figure  4.1a)  with  that  of 
the  three-dimensional  plot,  Figure  4.1c,  we  see  that  beam  pattern  of  the  two- 
dimensional  graph  is  identical  in  the  0  =  0°  plane  of  the  three-dimensional  plot. 


Figure  4.1a  2-D  plot  from  the  computer  output. 

A  two-dimensional  plot  of  equation  (4.1),  of  F  values  between  -10.9  dB  and  -10.0 
dB  results  in  Figure  4. Id.  For  F  values  less  than  -90.0  dB  we  get  Figure  4.1e. 

In  the  case  where  F  <  -90  dB,  H  <  lO'^O  20  which  is  small  enough  to  be  assumed  to 
be  0,  |Cj|  ~  |c2|  or  0  =  nrt  n  =  0,1,2 — ,  and  equation  (2.40)  can  be  tan<J>  =  1  or  -1. 
From  these  resulting  equations  we  can  extract  O  =  45°  +  n7t,'2,  n  =  1,2  — .  As  seen 


in  Figure  4.1e,  when  F  <  -90  dB  the  asymptote  matches  the  nulls  in  the  0-<I> 
relationship  of  Figure  4.1c.  Note  the  nulls  at  0  -  0°,  180°  for  all  <D,  and  <D  «  45°. 
These  figures  of  constant  dB  level  aid  in  the  interpretation  of  Figure  4.1c.  The 
modified  P3D  program  (TRACE)  is  included  in  Appendix  F. 
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K=10.S5 

Dl=  0.5  D2=  05 
D3-  0.5  D4-  0.5 
VIEW  ANGLE—20.-60 
P2-  0  P3-18C  P4-180 
A=1  3=1  C=l. 


TRACK  (-10.  DB) 
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2  =  -90  DB 
K='.0  55 
Dl-05  D2-05 
D3-C5  34-0* 

P2=  C  F3=ieC  F4-10O 
A=1  B=1  C  =  1 


B.  CASE  2 

In  this  case,  a  triplet  array  is  investigated.  Referring  to  Figure  2.3,  the  fourth 
element  is  turned  ofT  and  its  opposing  element  is  moved  toward  the  center  of  the  array 
so  its  distances  is  half  that  of  the  two  elements.  Also,  the  source  2  is  180°  out  of  phase 
with  the  other  two  elements. 

When  kDj  -  kD;  -  5.2758. 

kD3  -  2.63"9,  kD4  -  0. 

<p,  -  Jt.  <p3  -  <p4  -  0. 

A  -  B  -  1.  C  -  0 
the  equation  (2.3")  can  be  simplified  as 


H(0.<D>  -  -  ; [cos< 2.6379  sinflcosO »):  +  [2  sin(5.2"58  sinOsinO) 

-  sin( 2.63"*9  sinBcos^)]*1  :a. 


If  values  are  calculated  from  the  equation  <4  2)  as  in  case  I.  then  comparisons  can  be 
made  between  these  calculated  values  and  the  two-dimensional  plot  as  seen  in  Figure 
4.2a  i for  0  m  0°.  20°  and  90°.  and  <P  *  90°,  we  get  F  ■  -9.5424.  -2. "431  and  -3.0  o 
respectively  from  hand-computation./.  Additionally,  the  two-dimensionai  plot 
corresponds  well  in  shape  with  Figure  4.2b  from  the  experiment  t  Ignore  the  ordering  of 
the  relative  phases  in  this  figures;  it  is  the  consequence  of  a  different  laDeiing 
convention). 

If  we  compare  the  shape  of  a  two-dimensional  plot  ( Figure  4.2a)  w-ith  that  of  the 
three-dimensionai  piot  (Figure  4.2o,  we  see  that  beam  pattern  ot  the  two-dimensionai 
graph  is  identical  in  the  <l>  *  90°  plane  of  the  three-dimensional  plot. 

For  a  contour  piot  of  this  case,  a  two-dimensical  piot  ol  equation  '  -4  2 1  ot  F 
values  between  -5.9  dB  and  -5.0  dB  results  Figure  4. 2d.  For  F  values  between  -10.9  dB 
and  -10.0  dB,  we  get  Figure  4.2e  For  F  values  less  than  -35  0  dB.  we  get  Figure  4 . 2 f . 
From  equation  '4.2).  we  get  F  =  -35.8~0  for  0  =  143"  and  =  10°  from  hand- 
calculation.  This  point  is  seen  in  Figure  4  2f.  Successive  examination  of  Figure  4, 2d 
through  4.2f  shows  how  the  contour  lines  migrate  as  lower  and  lower  values  of  F  are 
examined  These  figures  of  constant  dB  level  aid  in  the  interpretation  of  Figure  4.2c. 
The  modified  P3D  program  (TRACE)  is  included  in  Appendix  F 


K=l?586 
Dl=  0.3  D2=  0.3 
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TRACK  (-5.0  DB) 
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C.  CASE  3 

In  this  case,  a  doublet  array  is  investigated.  Referring  to  Figure  4.2a,  this  is  the 
same  as  the  first  case  except  that  the  third  and  fourth  elements  are  turned  off  and  all 
elements  are  in  phase. 

For  kDj  =  kD,  =  kD3  =  kD4  =  5.275, 

<P2  =  <P3  -  <P4  =  0, 

A-l,  B-C-0 
equation  (2.37)  can  be  simplified  as 

H(0,O)  *  |cos(5.275  sin0sin<I>)| 


As  in  cases  1  and  2,  for  0  =  0°,  35°,  90°  and  <P  =  90°,  we  get  F  *  0,  -0.058  and 
-5.459  respectively  from  hand-computation.  These  points  are  found  in  Figure  4.3a  and, 
the  two-dimensional  plot  again  corresponds  well  in  shape  with  Figure  4.3b  from  the 
experiment. 

If  we  compare  the  shape  of  a  two-dimensional  plot  (Figure  4.3a)  with  that  of  the 
three-dimensional  plot  (Figure  4.3c  below),  we  see  that  the  beam  pattern  of  the  two- 
dimensional  graph  is  identical  in  the  <D  =  0°  plane  of  the  three-dimensional  plot. 

For  a  constant  value  of  H(0,<1>),  we  can  follow  the  same  procedure  to  get  the 
equation  for  contour  plots  as  in  case  1  and  2.  The  resultant  figures  are  Figures  4.3d, 
4.3e  and  4.3f.  Figure  4.3d  is  drawn  for  the  F  values  between  -5.9  dB  and  -5.0  dB  and 
Figure  4.3e  is  for  the  F  values  between  -8.9  dB  and  -8.0  dB.  Figure  4.3f  is  for  F  values 
less  than  -35.0  uB. 
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Figure  4.3a  2-D  plot  from  the  computer  output 


Figure  4.3b  2-D  plot  from  the  experiment 
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TRACK  (-  5.0  DB) 


Figure  4.3d  Track  (-5  dB). 


TRACK  (-35.0  DB) 


THETAOEC) 


2  =  -360  DB 
K=:0.55 
Dl-0  5  D2-C5 
D3«:5  34=)  f 
F2=  0  ?3*=  )  P4= 
A=1  B=0  0=0 


Figure  4. 3f  Track  (-35  d B ). 
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V.  EVALUATION 


In  the  last  chapter,  we  looked  at  three  simple  cases.  These  established  the 
accuracy  of  the  five  computer  programs.  For  example,  in  comparisons  between  hand 
calculation  and  computer  output,  both  sets  of  values  were  proven  nearly  identical. 
Again,  due  to  the  complexity  of  calculations,  not  all  values  were  hand-calculated. 

When  a  two-dimensional  circular  graph  was  compared  with  the  selected  hand- 
computed  values,  both  groups  of  values  were  again  virtually  identical.  As  another  test 
of  the  validity  of  the  five  computer  programs,  some  of  the  computer  generated  two- 
dimensional  plots  were  compared  against  those  from  an  experiment  (Ref.  9:  p.18). 
Visual  comparisons  of  the  two  types  of  plots  were  similar  in  shape,  but  values  in  dB 
were  somewhat  different.  These  differences  might  be  caused  by  the  normalization  of 
the  angular  dependance  of  the  sound  pressure.  Also,  because  of  the  directivity  of  the 
real  sources,  there  is  a  slight  decrease  in  dB  level  of  the  experimental  plots  when  the 
array  pointed  away  from  the  receiver.  However,  they  are  small  enough  that  the 
comparisons  stiil  seem  quite  reasonable. 

Next,  two-dimensional  versus  three-dimensional  graphs  were  compared. 
Agreement  here  shows  the  programs  are  working  properly.  By  varying  the  range  of  0 
and  <J>  that  are  plotted,  and  also  by  vanng  the  view  angles,  values  and  the  shape  of  the 
curve  can  be  visualized  without  great  difficulty. 

Tracks  of  certain  directivity  values  were  drawn  using  the  modified  P3D  program 
(TRACE)  and  were  compared  to  the  three-dimensional  plots.  The  shape  of  the 
resuiting  tracks  corresponded  well  with  those  of  the  three-dimensional  patterns,  again 
showing  proper  functioning  of  the  programs. 

With  all  the  comparisons  made  in  Chapter  IV,  the  five  computer  programs 
developed  in  this  thesis  were  proven  effective  in  predicting  the  directivity  pattern  in  all 
directions  of  space  with  or  without  differences  in  the  strength  of  the  sources,  inter- 
source  distances  and  phases. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

The  following  conclusions  are  offered. 

•  The  P2DICIR,  P2DPCIR,  P2DIFIX.  P2DPF1X  programs  successfully 
predict  the  directivity  beam  pattern  in  two  dimensions. 

•  The  P3D  program  predicts  and  visualizes  the  directivity  beam  pattern  in  the 
three-dimension  space. 

•  They  allow  the  presentation  of  complete  beam  patterns,  which  contain 
considerable  information,  in  a  succinct,  easy  to  comprehend  and  interpret  form. 

RECOMMENDATIONS 

The  following  recommendations  and  areas  of  further  investigation  is  suggested. 

•  Further  studies  for  computer  programs  to  visualize  the  beam  pattern  in  the 
x-y-z  space  are  needed  to  eliminate  hidden  lines  in  three-dimensional  plotting. 

•  Transcendental  equations  need  to  be  solved  to  get  the  contour  plot  for  fixed 
F  values. 


APPENDIX  A 

BASIC  PROGRAMS  (HP  COMPUTER) 


THE 

OBJECTIVE  OF  THESE  PROGRAMS  IS  TO  CALCULATE  ANO 

ORAM  THE  DIRECTIVITY  PATTERN  IN  TWO-DIMENSION. 

THIS  IS  WRITTEN  IN  BASIC  LANGUAGE  ANO  PLOTTER  IS  HP7090A. 

2  : 

OIRECTIVITY  IN  dB 

P2 

PHASE  1  DEGREE  )  OF  SOURCE  A 

PS 

PHASE  l DEGREE  )  OF  SOURCE  B 

P4 

PHASE  ( DEGREE )  OF  SOURCE  C 

P5 

VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN 

P6 

RADIAL  VALUE  OF  PI 

P7 

RADIAL  VALUE  OF  P2 

P8 

RADIAL  VALUE  OF  P3 

K  : 

WAVE  NUMBER 

D  : 

INTER-SOURCE  DISTANCE 

P  : 

PHI  IN  DEGREE 

I  : 

THETA  IN  OEGREE 

PI 

3.141S917 

HP  PROGRAM  FOR  THE  RECTANGULAR  GRAPH 

INPUT  PARAMETERS 


5  OPTION  BASE 

6  RAO 

10  DIM  Z\  360  ) 


20 

OISP  1 

•WHAT 

IS 

THE 

VALUE 

OF 

KO?  KO*  "  ! 

K0-K«0 

30 

INPUT 

KD 

60 

OISP  • 

■WHAT 

IS 

THE 

VALUE 

OF 

P2(  OEG )?  P2* 

P2*PHASE 

OF 

A 

61 

INPUT 

P2 

62 

OISP  ■ 

•WHAT 

IS 

THE 

VALUE 

OF 

P3(  DEG )?  P3* 

•*  1 

PS^PHASE 

OF 

B 

63 

INPUT 

P3 

64 

OISP  ' 

•WHAT 

IS 

THE 

VALUE 

OF 

P3I0EGI?  P4® 

"  ! 

P4*PHASE 

OF 

C 

65 

INPUT 

P4 

66 

DISP  1 

■WHAT 

IS 

THE 

VALUE 

OF 

A?  A*  ’• 

67 

INPUT 

A 

68 

OISP  1 

•WHAT 

IS 

THE 

VALUE 

OF 

B?  B=  '* 

69 

INPUT 

B 

70 

OISP  • 

•WHAT 

IS 

THE 

VALUE 

OF 

C7  C=  ” 

71 

INPUT 

C 

72 

OISP  • 

•IS  P 

FIXED? 

(YES)ANsl  ,  NO | AN=0 ) 

AN*  M 

73 

INPUT 

AN 

74 

IF  AN= 

=1  THEN  DISP 

"WHAT 

IS 

THE  VALUE  OF 

P(DEG)? 

P*  "ELSE  GOTO 

76 

75 

INPUT 

P 

76 

DISP  1 

•IS  I 

FIXED? 

(YESlANSs 

* 1  ,  NO lANS=0 ) 

ANS* 

tl 

77 

INPUT 

ANS 

78  IF  ANS=1  THEN  OISP 

79  INPUT  I 


•HHAT  IS  THE  VALUE  OF  I(DES)?  1=  "ELSE  SOTO  82 


82  OISP 


83  IF  AN*0  ANO  ANS=0  THEN  SOTO  86 

84  IF  AN=0  AND  ANS*1  THEN  SOTO  90 

85  IF  AN*1  AND  ANS=0  THEN  GOTO  93 

86  OISP  »"K*0*"  »KD »"  ••,"*•• 


87  OISP  ”*PHASE  OF  A( OEG )«" >P2>"PHASE  OF  B( DEG )*" iP3 ("PHASE  OF  C( OEG )*" >P4,' 

88  OISP  >“A=" >A  >"  S3" »B»"C3" jC  >"»"  3  GOTO  99 


90  OISP 

91  DISP 

92  DISP 

93  DISP 


•*"»"K*0*"(K0, 


I(  OEG  )*"  »I 


’•PHASE  OF  AtOEG)3">P2j"PHASE  OF  Bl DEG  )=” »P3 »"PHASE  OF  Cl  DEG  )«"  »P4 


’ » "A*‘* » A , "  B3‘ 
1  »"K*03"  IKO  »* 


>B,"C3">C, 
I(  OEG  )*"  )I 


3  GOTO  103 


94  OISP  "*PHASE  OF  A( OEG  )*" (P2  > "PHASE  OF  B( DEG  )*" )P5> "PHASE  OF  C(DEG)>“>P4 

95  OISP  i"A*“ ( A ,"  B=” )B»"C=” }C  »"*”  3  GOTO  107 


100  OISP  ”  “ 

101  DISP  "  THETA” j” 

102  FOR  P>1  TO  360  3  GOTO 

PHI”,” 

no 

DIRECTIVITY” 

1U6 

104 

105 

DISP  "  ” 

OISP  ”  “v”  THETA” t" 

PHI",” 

DIRECTIVITY" 

106 

FOR  P*1  TO  360  3  GOTO 

115 

107 

108 
109 

DISP  “  ” 

DISP  ”  THETA”  >" 

PHI”,” 

DIRECTIVITY" 

MAIN  PROGRAM 


UO  FOR  1  =  1  TO  360 

ns  p5=pi  /iao 

116  P6*P2*P5  3  P7=P3*P5  3  P8*P4*P5 

117  9*P-90 

118  Q1=SIN  (0»P5) 

119  Il-SIN  ( I*P5  ) 

120  Q2=SIN  (P*PS) 

121  Hl=KD*Il*02/2 
121  H2=KD*Il*91/2 

127  W=COS  <H1!+A*C0S  < N1*P6  )+B*COS  <  W2-P7  )*C*COS  (N2+P8) 

128  G=SIN  ( Ml )-A*SIN  ( H1*P6  )*B*SIN  ( W2-P7  )-C*SIN  IW2+P8) 
130  H*SQR  ( H**2+G**2 )/< 1+A*B*C  ) 

140  Z(  P  )=20*LGT  (H) 

143  T=20*LGT  l  ABS  (SIN  <.5*KD*Iim 
150  DISP  USING  153  »  I, P, ZIP  I 

153  IMAGE  5X,3D . 0 ,5X ,30 . D , 15X ,S000 . ODD  > 15X ,SODD. DDD 

155  !  IF  Z(  I  X-20  THEN  1*1*.  1  3  GOTO  119 

156  •  I*  IP<  I  ) 

157  !  IF  ZtP)<-20  THEN  P=P+.l  3  GOTO  117 

158  !  P«  IP( P  >> 

IbO  IF  AN=1  AND  ANS=0  THEN  GOTO  163  ELSE  GOTO  161 

161  IF  AN=0  AND  ANS=1  THEN  GOTO  164  ELSE  GOTO  163 

162  NEXT  P  3  GOTO  96 

163  *  NEXT  I 

164  NEXT  P 


PLOTTING  WITH  HP7090A 


170  PLOTTER  IS  705 
180  CLEAR 
190  GCLEAR 
200  GRAME 

210  LOCATE  15,145,10,80 

220  SCALE  0,390,-240,0 

230  CLIP  0,360,0,-240 

240  AXES  5,5.0,-240,10,10,2 

250  LINE  TYPE  1 

260  GRID  20,-20 

270  FRAME 

280  LORG  6 

290  CSIZE  4 

300  FOR  1=0  TO  360  STEP  40 
310  MOVE  1,-241 

320  LABEL  I 

330  NEXT  I 
340  LORG  8 

350  FOR  I  =  -240  TO  0  STEP  40 

360  MOVE  -1,1 

370  LABEL  I 

380  NEXT  I 

390  CSIZE  5 

400  LORG  5 

410  MOVE  180,40 

420  CSIZE  * 

430  PEN  2 

440  LABEL  "PHI  VS  DIRECTIVITY” 

450  LORG  5 

460  MOVE  180,20 

470  CSIZE  4 

480  LABEL  "(THETA  *  45  DEGREE)" 

490  CSIZE  5 

500  LORG  5 

510  MOVE  180,-26- 

520  CSIZE  4 


, 


530  LABEL  "PHI (DEG)" 

540  LOIR  90 
550  HOVE  -30,-120 
560  DEG 
570  LDIR  90 
580  CSIZE  4 

590  LABEL  "DIRECTIVITY* dB )" 

600  CSIZE  5 

610  MOVE  0,0 

620  FOR  P*1  TO  360 

630  DRAM  P,Z(P) 

640  NEXT  P 
650  PEN  0 
700  END 
117392 


INPUT  PARAMETERS 


5  OPTION  BASE  1 

6  RAD 

10  DIM  XI 360) 

11  DIM  Y<  360  ) 

12  DIM  Z( 360  ) 

20  DISP  "WHAT  IS  THE  VALUE  OF  KD?  KD*  "  !  KD=K*D 

30  INPUT  KD 

60  DISP  "WHAT  IS  THE  VALUE  OF  P2IDEG)?  P2*  "!  P2=PHASE  OF  A 

61  INPUT  P2 

62  DISP  "WHAT  IS  THE  VALUE  OF  P3IDEG)?  P3=  "  !  P3=PHASE  OF  B 

63  INPUT  P3 

64  DISP  "WHAT  IS  THE  VALUE  OF  P4I0EG)?  P 4*  "  •  P4=PHASE  OF  C 

65  INPUT  P4 

66  DISP  "WHAT  IS  THE  VALUE  OF  A?  A*  " 

67  INPUT  A 

68  DISP  "WHAT  IS  THE  VALUE  OF  B?  B*  “ 

69  INPUT  B 

70  DISP  "WHAT  IS  THE  VALUE  OF  C?  C*  " 

71  INPUT  C 

72  DISP  "IS  P  FIXED?  !YESvAN*l  ,  N0>AN*0)  AN*  " 

73  INPUT  AN 

74  IF  AN*1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  P( DEG )?  P*  "ELSE  GOTO  76 

75  INPUT  P 

76  DISP  "IS  I  FIXED?  <YES>ANS*1  ,  NO|ANS=0 )  ANS*  " 

77  INPUT  ANS 

78  IF  ANS=1  THEN  DISP  "WHAT  IS  THE  VALUE  OF  II  DEG)?  I*  "ELSE  GOTO  82 

79  INPUT  I 

82  DISP  "IHHHHHUHHHHHHHHHHHHHHHHHHHHHHHHUHHHHHHHHHUHHHHHWHHHHHt" 

83  IF  AN=0  AND  ANS=0  THEN  GOTO  86 

84  IF  AN=0  AND  ANS=1  THEN  GOTO  90 

85  IF  AN=1  ANO  ANS=0  THEN  GOTO  93 

86  DISP  "*"l"K*0="»KD," 

87  DISP  "*PHASE  OF  A I  DEG  Is" »P2 ) "PHASE  OF  B1  DEG  )*" >P3 > "PHASE  OF  Cl  DEG  l=" >P4,"*" 

88  DISP  i "A=" j A , "  B=" |8  «"C  =  " jC S  GOTO  99 

90  DISP  >"K*D=" >KD ,"  II  DEG  )*" »I ," 

91  DISP  "*PHASE  OF  Al DEG  )*"  \P2 ("PHASE  OF  B<  DEG  )*" jP3  > "PHASE  OF  Cl OEG  )*" >P4,"»" 

92  OISP  >"A=" i A,"  B=" >B,"C*" >C S  GOTO  103 

93  DISP  "»"("K*0="jK0,"  II  DEG  )  =  " »I ," 

94  DISP  "*PHASE  OF  Al DEG  )=" >PZ ("PHASE  OF  B( DEG  )*" >P3  J "PHASE  OF  Cl  DEG  )=" jP4,"*" 

95  DISP  "»" >"A=" jA,"  B="jB,"C=">C,"»"  S  GOTO  107 

99  OISP  "xhhhhhhhhhhhhhhhhhhhhhhhhhhhuhhuhhhhhhhhhhhhhhhhhhhhi" 

100  OISP  "  " 

101  DISP  "  "J"  THETA" j"  PHI","  DIRECTIVITY" 

102  FOR  P=1  TO  360  3  GOTO  110 

’.03  OISP  . . . . . . . . 

104  DISP  "  " 

105  DISP  "  ">"  THETA" >"  PHI","  DIRECTIVITY" 

106  FOR  P=1  TO  360  3  GOTO  115 

107  DISP  "»»««« UK  . . . 

108  DISP  "  " 

109  OISP  "  " >"  THETA" J"  PHI","  DIRECTIVITY" 


MAIN  PROGRAM 


110  FOR  1*1  TO  360 
115  P5*PI  /180 


53 


3 


116  P6=P2*P5  3  P7=P3*P5  3  P8=P4«P5 

117  Q*P-90 

118  GlsSIN  (9»P5) 

119  I1=SIN  (I*P5) 

120  Q2=SIN  (P*P5) 

121  Wl=KD*Il*<32/2 
121  W2=KD*Il*Gl/2 

127  H=COS  <W1)+A«C0S  ( H1+P6 )+B*COS  ( W2-P7  )+C*COS  (W2+P8) 

128  G=SIN  (H1)-A*SIN  ( H1+P6  )*B«SIN  ( W2-P7  )-C*SIN  (W2«-P8) 
130  H=SQR  ( W*h»2+G#*2 )/( 1+A+B+C ) 

140  Z(I)=20*LGT  (HJ  +  240 

141  X(I  )=Z(I  )*COS  ( I*P5  ) 

142  Y(  I  )=Z(  I  )*SIN  <I«PS) 

143  T=20*LGT  ( ABS  (SIN  <.5»K0*I1))) 

ISO  OISP  USING  1S3  i  I,P,Z(P) 

153  IMAGE  5X  >SODD . ODD  » 15X.SDDD . 000 

155  !  IF  Z(  IK-20  THEN  I=I+.l  3  GOTO  119 

156  !  1=  IP(I) 

157  •  IF  ZlPK-20  THEN  P=P+.l  3  GOTO  117 

158  !  P=  IPtP)) 

160  IF  AN=1  AND  ANS=0  THEN  GOTO  163  ELSE  GOTO  161 

161  IF  AN=0  AND  ANS=1  THEN  GOTO  164  ELSE  GOTO  163 

162  NEXT  P  3  GOTO  96 

163  NEXT  I 

164  !  NEXT  P 


PLOTTING  WITH  HP7090A 


APPENDIX  B 

P2DIFIX,  P2DPFIX,  P2DICIR  AND  P2DPCIR 


WHHHHHHHHHHHHHHHHHHHHHt  INFORMATION  IHHHHHHHHHHHHHHHHHHHHHHHHt 
* 

*  THE  OBJECTIVE  OF  THESE  PROGRAMS  IS  TO  CALCULATE  AND 

m  DRAW  THE  DIRECTIVITY  PATTERN  IN  TWO-DIMENSION. 

*  THE  DRAWING  METHOD  IS  THE  SUBROUTINE  CURVE  OF  DISSPLA. 

«  1(0)  AND  P(<P)  CAN  BE  USED  INTERCHANGEABLY  FOR 

*  P20PFIX  AND  P2DICIR  PROGRAMS. 

« 

IHHHHHHHHHHHHHHHHHHHHHt  VARIABLE  DEFINITION  IHHHHHHHHUHHHHHHHHt 


* 

Z  : 

DIRECTIVITY  IN  dB 

« 

PI 

PHASE  (DEGREE)  OF  SOURCE  A 

* 

P2 

PHASE  (DEGREE)  OF  SOURCE  B 

* 

P3 

PHASE  (DEGREE)  OF  SOURCE  C 

tt 

P4 

VARIABLE  TO  CHANGE  DEGREE  TO 

RADIAN 

* 

PS 

RADIAL  VALUE  OF  PI 

« 

P6 

RADIAL  VALUE  OF  P2 

* 

P7 

RADIAL  VALUE  OF  P3 

* 

K  : 

HAVE  NUMBER 

« 

D1 

DISTANCE  FROM  THE  GEOMETRICAL 

CENTER 

TO 

THE 

SOURCE 

1 

* 

D2 

DISTANCE  FROM  THE  GEOMETRICAL 

CENTER 

TO 

THE 

SOURCE 

A 

« 

03 

DISTANCE  FROM  THE  GEOMETRICAL 

CENTER 

TO 

THE 

SOURCE 

B 

* 

D4 

DISTANCE  FROM  THE  GEOMETRICAL 

CENTER 

TO 

THE 

SOURCE 

C 

* 

P  : 

PHI  IN  DEGREE 

* 

I  : 

THETA  IN  DEGREE 

* 

PH 

PHI  IN  RADIAN 

* 

TH 

THETA  IN  RADIAN 

• 

PHI 

:  DUMMY  VARIABLE  OF  PHI 

*  THETA  :  DUMMY  VARIABLE  OF  THETA 

«  PI  :  3.1415917 

* 

tWHHHHHHHHHHHHHHHHHHHHt  VARIABLE  DECLARATION  WHHHHHMHHHHHHHH 
C 

c _ _ _ _ 

PROGRAM  FOR  THE  TWO-DIMENSIONAL  RECTANGULAR  GRAPH 


C 

DIMENSION  PHI ( 361 ) , THETA! 361 )  ,Z( 360  ) 

REAL  P4, P5, P6,P7,H1,H2,W3,H4,W,G,H,K, D1,D2,D3,D4, PH, TH, PI, a, 6, c 
INTEGER  P,P1,P2,P3,I 
C 

c - 

C  INPUT  PARAMETERS 

C - 

c 

PI=3. 1415917 

K=5.5 

Dl*l. 

02=1. 

D3=l. 

04*1. 

PI*  180 
P2*  0 
P3*  180 
A  *  1. 

B  *  1. 

C  *  1. 

I  *  25:  USE  P  FOR  P2DPFIX 
C 

C - 

C  MAIN  PROGRAM 

C - 

C 

P4*PI/180 
P5=P1*P4 
P6*P2*P4 
P7*P3*P4 

DO  100  1*1,360:  MOVE  C  TO  LINE  NO.  76  FOR  P2DPFIX 
DO  200  P*l,360 


c 


1 

zoo 

100 

c 

c — 
c 

c — 
c 


c 


c 

c 


TH*I*P4 

PH=P«P4 

H1>K«01«SXN( PH )*SIN( TH  ) 

H2«K»02*SIN( PH )*SIN( TH )_P5 
K3*K«D3*SIN< TH)*COS( PH )-P6 
H4=K*D4*SIN( TH  )#COS( PH  )+P7 
H*COS( HI )+A*COS< HZ )+B*COS( H3 )+C*COS( H4 ) 
G=SIN(  HI  )-A*SIN(  H2  )-B*SIN(  H3  HC#SIN(  H4 ) 
H=SQRT( H*H+G*G )/( 1+A+B+C ) 

Z( P  )=20»ALOG10< H ) 

IF(Z(P).LE.-40.  )  THEN 
Z( P )*-40 

ELSE 

Z(P)»Z(P) 

END  IF 
THETA!  P)=I 
PHI! P  )*P 

HRITE!  8»1 )  PHI!  P), ZIP) 

FORMAT! 2X »2( E10 . 3  >SX  ) ) 

CONTINUE 

CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK616 


CALL  MEOBUF 
CALL  COMPRS 

CALL  TEK618:  MOVE  C  TO  DRAH  USING  TEK618 
CALL  PHYSORI 1. ,1.1 
CALL  PAGE! 15. 5, 1Z.  ) 

CALL  AREAZO! 10.5,9.1) 

CALL  FRAME 
CALL  COMP LX 
CALL  HEIGHT (.2) 

CALL  XNAME! 'PHI (DEG)  *',10) 

CALL  YNAME!  'DIRECTIVITY!  OB  )*  ‘  ,15) 

CALL  YTICKS! 8 ) 

CALL  XTICKSI 8 ) 

CALL  GRAF! 0 ,40,360,-40,4,0 ) 

CALL  GRID! 2,2 ) 

CALL  CURVE! PHI ,Z, 359,0) 

CALL  MESSAG! 'THETA=25  DEGREE ' ,15,11. ,8.5 ) 

CALL  MESSAG! 'K=5. S' ,5,11. ,8.2) 

CALL  MESSAG! 'Dl*l.  DZ=1. ' ,11 ,11 . ,7. 9) 

CALL  MESSAG! '03=1.  D4«l. • ,11 ,11. ,7. 6 ) 

CALL  MESSAG! 'P2=180  P3=0  P4*180 ' ,18,11 . ,7.3 ) 
CALL  MESSAG! 'A=l.  B=l.  C=l. ' ,14,11. ,7.0 ) 

CALL  RESET! 'HEIGHT' ) 

CALL  HEIGHT! .4) 

CALL  MESSAG! 'PHI  VS.  DIRECTIVITY' ,19,3.1,9.8) 
CALL  RESET! 'ALL'  ) 

CALL  NOCHEK 
CALL  ENDGR(O) 

CALL  PHYSORI 11.85,7.6) 

CALL  AREAZD! 3. 5,2.5 ) 

CALL  FRAME 
CALL  ENDPL(O) 

STOP 

END 


PROGRAM  FOR  THE  CIRCULAR  THO-DIMENSIONAL  GRAPH 


C 

DIMENSION  PHI! 360  )  ,R( 360 ) ,THETAf  360  )  ,Z( 360 )  ,X( 4 ) ,Y<  4 ) 

REAL  P4 ,P5,P6 ,P7 ,H1 ,HZ ,H3 ,H4,H,G ,H,K ,D1 ,02 ,D3 ,D4 ,PH,TH,PI ,a ,b,c 
INTEGER  P  »P1  ,P2  ,P3  >1 >Q 
C 

C - 

C 

C - 

C 

PI*3. 1415917 
K=5.5 
01=1. 

02=1. 

D3  =  l. 

D4=l. 

Pl=  180 
PZ=  0 
P3=  180 


INPUT  PARAMETERS 


ooooo 


B  «  1. 

C  ■  1. 

I  »  25:  USE  P  FOR  P20PCZR 


MAIN  PROGRAM 


P4*PX/180 
P5*P1»P4 
P6=P2*P4 
P7*P3«P4 

C  DO  100  I >0,360:  MOVE  C  TO  LINE  NO.  174  FOR  P2DICIR 
DO  200  P*0 ,360 
TH*I»P4 
PH*P»P4 

H1=K*D1«SIN( PH )*COS(TH ) 

W2*K*02»SIN!  PH  )*COS<  TH  )+P5 
H3=K*D3«SIN( TH )*COS< PH )>P6 
H4=K*D4«SIN(  TH  )«C0S(  PH  )+P7 
W=COS(  HI  )+a*COSIH2  )+B*COS!W3  )+C*COS<  W4  ) 

G=SIN( HI  )-A*SIN( H2  )+B«SIN( W3 )-C*SIN( H4 ) 

H=SGRT<  W*W+G«G )/! 1+A+B+C ) 

Z(P)*20«ALOG10(H) 

IF(Z( P  ) . LE . -40 . )  THEN 
Z(P)*0. 

ELSE 

Z(  P  )*40 .  ♦Zl  P ) 

END  IF 
THETA! P  )*I 
PHI !P)*P 

WRITE (8,1)  PHK  P ) ,Z( P  ) 

1  FORMAT!  2X,Z<  E10. 3, 5X)) 

200  CONTINUE 

100  CONTINUE 
C 

C - 

C  A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 

C - 

C 

CALL  MEOBUF 
CALL  COMPRS 

C  CALL  TEK618:  MOVE  C  TO  THE  UPPER  LINE  WHEN  DRAWING  WITH  TEK618 
CALL  PHYSOR! 1.5,1.  ) 

CALL  PAGE! 15.5,12.  ) 

CALL  AREA2DI8. ,8.  ) 

CALL  COMP LX 
CALL  HEIGHT!. 2) 

CALL  POLAR! 3.1415917/180. ,10. ,4. ,4. ) 

CALL  GRIO! 1,1) 

CALL  CURVE! PHI ,Z, 360,0) 

CALL  MESSAG! 'THETA*25  OEGREE 1 ,15, 9. 7,7. 7 ) 

CALL  MESSAG! ’K*5.5' ,5, 9. 7, 7. 4) 

CALL  MESSAG! ‘01*1.  02*1 . • ,11 , 9 . 7 ,7 . 1 ) 

CALL  MESSAG! ’D3=l.  D4=l . • ,11 ,9. 7,6 .8  ) 

CALL  MESSAG! 'P2*180  PS*0  P4=180 ' , 18 ,9. 7,6.5 ) 

CALL  MESSAG! ’A*l.  B*l.  C*l. ' ,14,9. 7,6 . 2 ) 

CALL  MESSAG! 'O', 1,3. 9, 8.1) 

CALL  MESSAG! ’-10’ ,3,3.7,7.11 
CALL  MESSAG!  ’-20’ ,3, 3. 7,6.1  ) 

CALL  MESSAG! ’-30' ,3, 3. 7, 5.1) 

CALL  MESSAG! 1 -40 1 ,3, 3. 7,4.1 ) 

CALL  MESSAG!  *0! DEG)', 6, 7. 6, 3. 9) 

CALL  MESSAG! ’180' ,3,0. ,3.9) 

CALL  RESET! ’HEIGHT'  ) 

CALL  HEIGHT!. 4) 

CALL  MESSAG! ’PHI  VS.  DIRECTIVITY ’ ,19,3. 0,9.5 ) 

CALL  RESET! ’ALL’  ) 

CALL  NOCHEK 
CALL  ENOGRIO) 

CALL  PHYSOR! 10.95,6.8) 

CALL  AREAZD! 3 . 5 ,2. 5 ) 

CALL  FRAME 
CALL  DONEPL 
STOP 
END 
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jl  _  A.  1  .La 


J* 


APPENDIX  C 

PROGRAM  FOR  THE  3-D  GRAPH  USING  SURMAT  ABOUT  X-Y  PLANE 


wiHUMWwwMimMmiK  INFORMATION  »«nm»nin<iiimmmmnHim<  *»  .hhhhhh 

THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAW 
THE  DIRECTIVITY  PATTERN  IN  THREE -DIMENSION. 

THE  DRAWING  METHOD  IS  THE  SUBROUTINE  SURMAT  OF  OISSPLA. 

•********»*<hhh>  VARIABLE  DEFINITION  MMMMimimiDnmMimimximm 

Z  :  DIRECTIVITY  IN  dB 

ZMM  :  DOUBLE  MATRIX  TO  GET  X-Y  PLANE 

ZM  :  DUMMY  VARIABLE  OF  Z 

PI  :  PHASE  (DEGREE)  OF  SOURCE  A 

P2  :  PHASE  (DEGREE)  OF  SOURCE  B 

P3  :  PHASE  ( DEGREE )  OF  SOURCE  C 

P4  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN 

PS  :  RADIAL  VALUE  OF  PI 

P6  :  RADIAL  VALUE  OF  P2 

P7  :  RAOIAL  VALUE  OF  P3 

K  :  NAVE  NUMBER 

01  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1 

02  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A 

03  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B 

D4  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C 

PHI  :  PHI  IN  DEGREE 
I  :  THETA  IN  DEGREE 
PH  :  PHI  IN  RADIAN 
TH  :  THETA  IN  RADIAN 
PI  :  3. 14 1591 7 

mum  ■  «ii  Minimum  VARIABLE  DECLARATION  mmmmmmmmm  mm  m  m  turn  mm  m 
DIMENSION  Z( 8100 ) ,ZMM( 90  >  90  ) 

REAL  P,T,R.X,Y,ZM,a,b,c 

REAL  P4,P5,P6  ,P7 »W1 ,W2  >W3 »N4 ,N,G,H,K,D1,02,D3,04,PH,TH,PI 
INTEGER  PHI,P1,P2,P3,2 


PROGRAM  FOR  THE  UPPER  HEMISPHERE 


C  A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 

c  . 

CALL  COMPRS 

C  CALL  TEK616:  MOVE  C  TO  THE  UPPER  LINE  WHEN  DRAWING  WITH  TEK618 
CALL  RESET  (3HALL) 

CALL  PAGE!  15. ,9. ) 

CALL  PHYSORI 1. ,3.  ) 

CALL  SCMPLX 

CALL  AREA2DI  14. ,8. ) 

CALL  X3NAMEI  1H  ,1 ) 

CALL  Y3NAMEI 1H  ,1) 

CALL  Z3NAME11H  ,1) 

CALL  VOLM3D  ( 17. ,17. ,17. ) 

CALL  VUANGLI +20 . , *20 . ,500 . ) 

CALL  INTAXS 
CALL  ZAXANG  (-90. ) 

CALL  GRAF3D( 0. ,10. ,80. ,0. ,10. ,80. ,0. ,10. ,40. ) 

CALL  BLSUR 

CALL  BGNMATI  90,90) 


C  INPUT  PARAMETERS 

C  —  —  —  —  - . — 

PI*3. 1415917 
MPI 


MAIN  PROGRAM 


P4*PI/180 
P5*P1*P4 
P4*P2*P4 
P7*P3*P4 
DO  200  I*-90  >0 
TH=I*P4 

00  100  PHI*1 *340 *4 
PH*PHI*P4 

M1*K*D1*SIN( PH  )*SIN( TH ) 

N2*K*D2*SIN(  PH  )*SIN(  TH  l-PS 
H3*K»D3*SIN!  TH  )«C0SI  PH  )-P4 
W4*K*D4*SIN<  TH  )*COS(  PH  )*P7 
NsCOSIHl  )*AMC0S(H2  )*B*COSI  H3  UC*C0S1  H4  ) 
G*SIN<  HI )-A*SIN( M2 )-B*SINl H3 )+C«SIN< H4  ) 
H*SQRTIW*H+G«G  )/( 1*A*B*C ) 

2(P HI  )*20*AL0C  10(H) 

IFIZlPHI  1.LE.-40.  )  THEN 
2IPHI  1*0. 

ELSE 

21  PHI )*Z( PHI  ) 

END  IF 

R*40 . +Z<  PHI  I 
P«  FLOAT! PHI  I 
T * FLOAT ( I  I 

X*R*COS<  TH  IPCOS!  PH  1*40  . 

Y*R«C0S1  TH  l»SINI  PH  l*<*0  . 

2M«-R»SINI TH  ■ 

HRITE  18.1  I  X.Y  .2* 

1  FORMAT!  2X.JI  £10  S.5X  I  i 

CALL  SCTHATiX.t  .ZM.1.0  ' 

100  CONTINUE 

200  CONTINUE 
C 


A  PROGRAM  TOP  PLOTTING  l<  SHEBPAOR  •<  ’fn.H 


CALL 

EWNATi  ZW.O 

CALL 

CALL 

SURMA' 

CASH 

T**  .  .  *0  . 

%<r 

CAL. 

CAL. 

HARFtP- 

COMPL* 

CAL. 

HEIGHT 

; 

CAL. 

SURvlS 

»C'  'TH 

CA.. 

HESSAG 

sjp»  : ' : 

-  i 

CAL. 

»E  SE  ’ 

C  OMP.  • 

CAL. 

•ESI’ 

he  ;  sh* 

CAl.  MNMA'  «  *■ 


DS«1. 

D4»l. 

PI*  X80 
P2*  0 
P3*  180 
A  «  1. 

8*0. 

C  ■  0. 

C 

C - 

C  MAIN  PROGRAM 

C - 

C 

P4*PI/180 

P5*P1*P4 

P6*P2»P4 

P7*P3*P4 

DO  250  I *-180, -90 
TH*I«P4 

DO  ISO  PHI*1,360,4 
PH*PHI*P4 

H1=K*D1«SIN< PH )«COS( TH ) 

H2=K*D2*SIN( PH )*COS(  TH )+P5 
H3=K*D3#SIN( TH )#COS<  PH )-P6 
H4*K*04»SIN( TH  J*COS(  PH )*P7 
N*COS(  HI  )+A*COS(  H2  )+B*COS<  W3  )+C*COS(  H4  ) 
G=SIN( HI }-A*SIN( H2 )+B*SINl H3 )-C*SIN(  H4 ) 
H=S£3RT(  N*H+G#G  )/( 14A+B+C  ) 

Z( PHI  )*20*AL0G10( H  ) 

IF( Z( PHI ) .  LE . -40 . )  THEN 
Z(  PHI )*0. 

ELSE 

Z!  PHI )=Z( PHI ) 

END  IF 

R*40 . +Z( PHI  1 
P*  FLOAT (PHI) 

T*FLOAT( I  ) 

X=R*COS< TH  )«COS( PH  )*40. 

Y=R*COS(  TH  )*SIN(  PH  )+40 . 

ZM=R«SIN( TH ) 

HRITE (8,1)  X,Y,ZM 
1  FORMAT! 2X,3(E10. 3, 5X ) ) 

CALL  GETMAT(X,Y,ZM,1,0) 

150  CONTINUE 

250  CONTINUE 

CALL  ENDMAT!  Z7H,0 ) 

CALL  SURMAT  (ZMM, 1,90, 1,90,0) 

CALL  DASH 
CALL  MARKER! 17 ) 

CALL  COMP LX 
CALL  HEIGHT! .2) 

CALL  SURVIS!  'BOTTOM'  ) 

CALL  ENDPL(O) 

CALL  DONE PL 
CHECKSUMS  1236  2 
STOP 


APPENDIX  D 

PROGRAM  FOR  THE  3-DIMENSIONAL  GRAPH  USING  VECTOR  DRAWING  METH 


hhhhhhhhhhhhhm  INFORMATION  «hhhhhhhhhhhhhhhhhhhhhhhhhhhhh» 

THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  ORAM 
THE  DIRECTIVITY  PATTERN  IN  THREE -DIMENSION. 

THE  DRAWING  METHOD  IS  THE  VECTOR  DRAWING  METHOD  OF  DISSPLA. 


. . .  VARIABLE  DEFINITION  . . . 

* 

Z  :  DIRECTIVITY  IN  dB  * 

XT  :  "TO"  POINT  OF  X  IN  DRAWING  VECTOR  * 

YT  :  "TO"  POINT  OF  Y  IN  DRAWING  VECTOR  * 

ZT  :  "TO”  POINT  OF  Z  IN  DRAWING  VECTOR  * 

XF  :  "FROM"  POINT  OF  X  IN  DRAWING  VECTOR  « 

YF  :  "FROM"  POINT  OF  Y  IN  DRAWING  VECTOR  * 

ZF  :  "FROM"  POINT  OF  Z  IN  DRAWING  VECTOR  * 

PI  :  PHASE  (DEGREE)  OF  SOURCE  A  * 

PZ  :  PHASE  (DEGREE)  OF  SOURCE  B  * 

PS  :  PHASE  (DEGREE)  OF  SOURCE  C  * 

P4  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  » 

PS  :  RADIAL  VALUE  OF  PI  • 

P6  :  RADIAL  VALUE  OF  P2  * 

P 7  :  RADIAL  VALUE  OF  PS  • 

K  :  NAVE  NUMBER  « 

01  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1  * 

DZ  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A  * 

DS  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B  * 

04  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C  * 

PHI  :  PHI  IN  DEGREE  • 

I  :  THETA  IN  DEGREE  » 

PH  :  PHI  IN  RADIAN  • 

TH  :  THETA  IN  RADIAN  » 

PI  :  3.1415917  * 


•  « 

VARIABLE  DECLARATION  »<hhhhhhhhhhhhhhhhhhhhhhhh> 
REAL  R,Z,XT,YT,XF,YF,ZF,ZT,»,b,e 

REAL  P4,P5,P6,P7,W1,W2,N3,M4,W,G,H,K,D1,D2,D3,D4,PH,TH,PI 
INTEGER  PHI,P1,P2,P3,I 
C 

DIMENSION  Z( 0 : 2500 ) ,ZMM( 46 ,46  >, VI 2000 ) ,U(  2000 ) 

REAL  P,T,R,X,Y,ZM,ZZM,A,B,C,V,U 

REAL  P4,P5,P6 »P7 »W1 ,W2 ,N3,W4,W,G,H,K,D1,D2,DS,D4,PH,TH,PI 
INTEGER  PHI ,P1,P2,PS,I,S 
C 

C - 

C  INPUT  PARAMETERS 

C - 

C 

PI*3. 1415917 
K=10 .55 
01*0.5 
D2*0 .5 
03*0.5 
D4*0.5 
PI*  0 
P2=180 
P3*180 
A  *  1. 

B  *  1. 

C  *  1. 

C 

C - 

C  MAIN  PROGRAM 

C - - - - - - - - - 

C 

P4*PI/180 

P5*P1*P4 

P6=PZ*P4 

P7*P3*P4 

S»1 
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*  *  *  *  * 


uuuou 


DO  200  1*0,90 
TH*I*P4 

DO  100  PHI *0, 180 
PH*PHI»P4 

H1*K*D1*SINI  PH  )*SIN(  TH ) 

H2*K*D2*SINI PH )*SIN( TH )-P 5 
H3*K*03*SIN(  TH  )«COS(  PH  )-P6 
H4*K*04*SIN(  TH  )«COS(  PH  )*P7 
H*COS(  HI  )+A*COS<  H2  )+B*COS<  H3  )*C*COSt  H4 ) 

G*SIN( HI  )-A*SIN< H2  )-B*SIN< H3 )+C»SIN<  H4 ) 

H*SQRT( H*H+G*G )/< 1+A+B+C ) 

Z( PHI  )=20»ALOG10( H  ) 

IF( ( Z( PHI )  .ST . -10. J . AND . ( Z( PHI ) . LT. -9.  ) )  THEN 
Z( PHI  )■  -30. 

ELSE 

Z(  PHI )«  Z(PHI) 

END  IF 
ZM*Z( PHI ) 

C  HRITEI 8>1 )  I.PHI.ZM 

IF  (ZM  . EQ.  -30.  )  THEN 
V(SJ  *  FLOAT  I  PHI) 

UIS)  *  FLOAT  (I  ) 

C  HRITEI 8,1 )  VIS), UIS), ZM 

S=S+1 
END  IF 

C  HRITEI 8,1)  I.PHI.Z 

Cl  FORMAT!  2X, 2(13, 5X),F10. 5) 

C  HRITEI 8,1)  AXI L ) ,BVI L ) ,ZM 

Cl  FORMAT I 2X,3(F10.3,5X)) 

100  CONTINUE 

200  CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPA  OR  BY  TEK618 


CALL  MEDBUF 
C  CALL  COMPRS 
CALL  TEK618 
CALL  PHYSORI 1 . ,1. ) 

CALL  PAGE! 15.5,12. ) 

CALL  AREA20I 10.5,9.1) 

CALL  FRAME 
CALL  COMP LX 
CALL  HEIGHT!. 2) 

CALL  XNAMEI 'PHIIDEG)  **,10) 

CALL  YNAMEI  ' THETA! DEG  )$ ' ,10  ) 

CALL  YTICKSI 2  ) 

CALL  XTICKSI  4 ) 

CALL  GRAF I  0,20,180,0,10,90) 

CALL  CURVE I V,U,S,-1 ) 

CALL  MESSAGI  *THETA*25  DEGREE ’ ,15 ,11 . ,8.5 ) 
CALL  MESSAGI *K=5. 5' ,5,11. ,8.2) 

CALL  MESSAGI '01=1.  02*1 . 1 ,11 ,11. ,7. 9 ) 

CALL  MESSAGI  1 03*1 .  D4*l .• ,11,11 . ,7.6  ) 

CALL  MESSAGI *P2*180  P3*0  P4*180* ,18,11. ,7.3  ) 
CALL  MESSAGI ’ A*1 .  B*l.  C*l. * ,14,11 . ,7. 0 ) 
CALL  RESET I  'HEIGHT1  ) 

CALL  HEIGHT! .4) 

CALL  MESSAGI  'TRACK! -20 . OB  ) 1 ,13  ,3 . 1 ,9.8  ) 

CALL  RESETI  'ALL'  ) 

CALL  NOCHEK 
CALL  ENDGRIO) 

CALL  PHYSORI 11.65,7.6) 

CALL  AREA2D! 3 .5,2 .5  ) 

CALL  FRAME 
CALL  ENOPLIO) 

STOP 


V  V  ' 


APPENDIX  E 
P3D 


THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AND  DRAM 
THE  DIRECTIVITY  PATTERN  IN  THREE -DIMENSION. 

THE  0 RAWING  METHOD  IS  THE  SUBROUTINE  SURMAT  OF  OISSPLA 
ABOUT  O-tP  PLANE. 

« 

* 

• 

* 

« 

* 

2  : 

DIRECTIVITY  IN  dB  „  _ 

» 

« 

Z791 

:  DOUBLE  MATRIX  TO  GET  U-<P  plan* 

* 

X  : 

DUMMY  VARIABLE  OF  THETA 

* 

Y  : 

DLfrPTY  VARIABLE  OF  PHI 

* 

ZM 

DUmY  VARIABLE  OF  DIRECTIVITY 

« 

PI 

PHASE  (DEGREE)  OF  SOURCE  A 

* 

P2 

PHASE  (DEGREE)  OF  SOURCE  B 

• 

PS 

PHASE  ( DEGREE  )  OF  SOURCE  C 

« 

P4 

VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN 

« 

PS 

RADIAL  VALUE  OF  PI 

» 

P6 

RADIAL  VALUE  OF  P2 

« 

P7 

RADIAL  VALUE  OF  P3 

* 

K  : 

HAVE  NUMBER 

• 

D1 

DISTANCE  FROM  THE  GEOMETRICAL  CENTER 

TO 

THE 

SOURCE 

1 

* 

D2 

DISTANCE  FROM  THE  GEOMETRICAL  CENTER 

TO 

THE 

SOURCE 

A 

• 

03 

DISTANCE  FROM  THE  GEOMETRICAL  CENTER 

TO 

THE 

SOURCE 

B 

* 

D4 

DISTANCE  FROM  THE  GEOMETRICAL  CENTER 

TO 

THE 

SOURCE 

C 

« 

PHI 

:  PHI  IN  DEGREE 

* 

I  : 

THETA  IN  DEGREE 

« 

PH 

PHI  IN  RADIAN 

* 

TH 

THETA  IN  RADIAN 

* 

PI 

3.1415917 

HHH 

* 

* 

DIMENSION  2( 2500  I  ,ZMMI 46  >46  ) 

REAL  X,Y,ZM,P4,P5,P6,P7,H1,W2,W3,H4,M,G,H,IC,D1,D2,D3,D4,PH,TH,PI 
ml  a,b,e 

INTEGER  PHI,P1,P2,P3,I 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 


CALL  COMPRS 

CALL  TEK618 :  MOVE  C  TO  THE  UPPER  LINE  HHEN  DRAWING  WITH  TEK618 
CALL  RESET  ( 3HALL  ) 

CALL  PHYSOR) 0.5,1. J 
CALL  PAGE) 15.5,12. I 
CALL  SCMPLX 
CALL  AREA20I 10.5,4.15) 

CALL  X3NAMEI  'THETAf DEGREE  )' ,13  ) 

CALL  Y3NAMEI *PHI( DEGREE >' ,11 > 

CALL  Z3NAMEI 1 OIRECTIVITY) DB ) ' ,15) 

CALL  VOLM30  1 17. ,17. ,17. ) 

CALL  VUANGLt *20. ,+20. ,500. ) 

CALL  INTAXS 
CALL  2AXANG  I -90. 1 

CALL  CRAF3DI0. ,5. ,90., 180. ,6., 360. ,-40. ,2. ,0. ) 

CALL  BLSUR 

CALL  BGNMATC 46,46) 


MAIN  PROGRAM 


.w  w  w-  V  v  w-j  -  w  w  w  t  »■  y  -w  ,  w  -  *. 


Pl«  180 
P2«  0 
P3*  180 
A  «  1. 

8  *  1. 

C  ■  1. 


P4*PI/180 
P5*P1*P4 
P6*P2»P4 
P7*P3*P4 
00  200  1*0,90,2 
TH*I*P4 

DO  100  PHI*180,360,4 
PH=PHI*P4 

W1*K»D1»SIN<  PH  l»SIN<  TH  ) 

W2*K«D2»SIN( PH  >*SIN< TH  )-P5 
W3*K*D3#SIN( TH  )*C0S( PH  )-P6 
W4=K*D4*SIN(  TH  >«C0S(  PH  )+P7 
W=COS(  HI  )+A*COS<  HZ  )+B*COS(  N3  )+C*COS(  M4  ) 
G*SIN( Ml  )-A*SIN( HZ  )-B*SIN( H3 )+C*SIN( M4  ) 
H*SQRTI N»H«-G*G  )/( 1*A*B*C  ) 

ZIPHI  )*20*ALOG10(H) 

IF(  Z(  PHI  1.LE.-40.  )  THEN 
Zt  PHI )*-40 . 

ELSE 

ZIPHI)*  ZIPHI) 

END  IF 

X*I 

Y*PHI 

ZM=Z< PHI  I 

WRITE!  8.1 )  X.Y.ZM 

FORMAT! ZX.3IE10. 3, 5X)) 

CALL  GETMAT(X,Y,ZM,1,0  I 
CONTINUE 
CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 


CHECKSUMS  12362 


ENDMATIZMM.O) 

SURMAT  <ZMM, 1,46, 1,46,0) 

DASH 

MARKER!  17) 

COMP LX 
HEIGHT! .2) 

SURVIS!  'BOTTOM'  ) 

MESSAG!  'K=5.5'  ,5,11.  ,8.5) 

MESSAG!  *01*1.  02*1. ' ,11,11. ,8.2) 

MESSAG! '03=1.  D4*l. ' ,11,11. ,7.9) 

MESSAG!  'VIEW  ANGLE*  20,  20  '  ,  18 , 11. ,7. 6  ) 
MESSAG! 'P2=180  P3=0  P4=180 ' ,18,11 . ,7. 3 ) 
MESSAG!  'A*l.  B*1 .  C*l. ' ,14,11. ,7. 0 ) 

RESET! 'HEIGHT' ) 

HEIGHT! .4 ) 

MESSAG! 'DIRECTIVITY  PATTERN' ,19 ,3 . 05,9.8 ) 
ENDGR! 0 ) 

PHYSOR! 11.35,7.6) 

AREA2D! 3. 5, 2. 5 ) 

FRAME 

OONEPL 


APPENDIX  F 
TRACE 


. .  INFORMATION  iHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHt 

• 

THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  AM}  DRAM  * 

THE  CONSTANT  dB  LEVEL  OF  THE  THREE-DIMENSIONAL  PLOT.  « 

THE  DRAWING  METHOD  IS  THE  SUBROUTINE  CURVE  OF  DISSPLA.  * 

• 

HHHHHHHHHHHHM  VARIABLE  DEFINITION  IHHHHHHHHHHHHHHHHHMHHHHHHHHt 

• 

Z  :  DIRECTIVITY  IN  e©  * 

v  :  ARRAY  TO  TRACK  U>  * 

u  :  ARRAY  TO  TRACK  V  * 

PI  :  PHASE  (DEGREE)  OF  SOURCE  A  • 

P2  :  PHASE  (DEGREE)  OF  SOURCE  B  * 

PS  :  PHASE  (DEGREE)  OF  SOURCE  C  • 

P4  :  VARIABLE  TO  CHANGE  DEGREE  TO  RADIAN  * 

PS  :  RADIAL  VALUE  OF  PX  • 

P6  :  RAOIAL  VALUE  OF  P2  » 

P7  :  RADIAL  VALUE  OF  PS  * 

K  ;  HAVE  NUMBER  * 

D1  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  1  * 

D2  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  A  « 

D3  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  B  * 

D4  :  DISTANCE  FROM  THE  GEOMETRICAL  CENTER  TO  THE  SOURCE  C  * 

PHI  :  PHI  IN  DEGREE  • 

I  :  THETA  IN  OEGREE  * 

PH  :  PHI  IN  RAOIAN  • 

TH  :  THETA  IN  RADIAN  * 

PI  :  S. 1415917  * 


HHHHHHHHHHHHHHHHI  VARIABLE  DECLARATION  «hhhmhhhhhhhhhhmh» 
DIMENSION  Z(  0 :  2500  )  ,Zht1(  46  ,46  ) ,  V(  2000  )  ,U(  2000  ) 

REAL  P,T,R,X,Y,ZM,ZZM,A,B,C,V,U 

REAL  P4,P5,P6,P7,W1,W2,W3,W4,W,G,H,K,D1,D2,DS,D4,PH,TH,PI 
INTEGER  PHI,P1,P2,P3,I,S 


INPUT  PARAMETERS 

W - . 

c 

PI*3. 1415917 
K*10 .55 
01*0.5 
02*0.5 
D3*0 .5 
04=0 . 5 
PI*  0 
P2*180 
P3=180 
A  *  1. 

B  *  1. 

C  «  1. 


MAIN  PROGRAM 


P4=PI/180 

P5*P1*P4 

P6*P2*P4 

P7*PS*P4 

S*1 

DO  200  1*0,90 
TH*I»P4 

DO  100  PHI *0,180 
PH=PHI*P4 

H1=K*D1*SIN(  PH  )*SIN(  TH  ) 

H2=K»D2*SINI PH  )*SIN( TH  )-P5 
HS*K*DS*SIN( TH  )*COS( PH  )-P6 
N4*K*D4«SIN( TH  )»COS( PH  )+P7 
H=COS(  HI  )+A*COS(  M2  )  ♦B*C0S1  H3  )+C*COS(  W4  ) 
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.  ' 


IT**1 


ooooo 


G*SIN( HI  )-A«SIN< HZ )-B*SIN< N3 )+C*SIN( H4 ) 
H=SQRT(H»W+G*G)/l  1+A+B+C  ) 

Zl PHI  )=20*AL0G10I H ) 

IF( ( Z( PHI )  . G7 . -10 . ) . ANO. ( Z( PHI ) . LT. -9. ) )  THEN 
Z( PHI  )*  -30. 

ELSE 

Z!PHI)«  Z(PHI) 

END  IF 
ZM=Z( PHI ) 

C  HRITE(8,1)  I,PHI,ZM 

IF  <ZM  . EQ.  -30.  )  THEN 
VIS)  *  FLOAT  (PHI ) 

U(S)  *  FLOAT  (I) 

C  WRITE (8,1 )  V( S )  ,U( S  )  ,ZM 

S*S+1 
END  IF 

C  WRITE ( 8,1 )  I ,PHI ,Z 

Cl  FORMAT(2X,2(I3,5X),F10.5) 

C  WRITE) 8,1 )  AX( L ) ,BY( L ) ,ZM 

Cl  FORMAT! 2X,3( F10.3,5X) ) 

100  CONTINUE 

200  CONTINUE 


A  PROGRAM  FOR  PLOTTING  BY  SHERPA  OR  BY  TEK618 


CALL  MEDBUF 

C  CALL  COMPRSlMOVE  C  TO  TEK619  TO  PRINT  USING  SHERPA 
CALL  TEK618 
CALL  PHYSORd.  ,1.  ) 

CALL  PAGE! 15.5,12.  I 
CALL  AREA2D! 10.5,9.1) 

CALL  FRAME 
CALL  COMP LX 
CALL  HEIGHT!. 2) 

CALL  XNAME!  'PHI! OEG)  *',10) 

CALL  YNAME!  'THETA! OEG  )* ' ,10  ) 

CALL  YTICKSI 2 ) 

CALL  XTICKS! 4 ) 

CALL  GRAF (0,20, 180 ,0,10,90) 

CALL  CURVE! V,U,S,-1) 

CALL  MESSAG! 'THETA=2S  DEGREE ' ,15,11. ,8.5 ) 

CALL  MESSAG! 'K*5. 5' ,5,11. ,8.2) 

CALL  MESSAG! 'Dl=l.  02=1. ' ,11 ,11. , 7. 9 ) 

CALL  MESSAG! '03=1.  04-1. ' ,11,11 . ,7.61 
CALL  MESSAG! 'P2=180  P3=0  P4=180 • ,18,11 . ,7.3 ) 

CALL  MESSAG!  'A=l.  8=1.  C=1 . ' ,14,11. ,7. 0  ) 

CALL  RESET!  'HEIGHT'  ) 

CALL  HEIGHT! .41 

CALL  MESSAG! 'TRACK! -20 .DB )' ,13,3 . 1,9.8 ) 

CALL  RESET!  'ALL'  ) 

CALL  NOCHEK 
CALL  ENDGR(O) 

CALL  PHYSOR! 11.85,7.6) 

CALL  AREA2D! 3. 5, 2.5 ) 

CALL  FRAME 
CALL  ENDPL(O) 

STOP 

END 
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APPENDIX  G 

LIST  OF  COMPUTER  OUTPUT 


•  SAMPLE  CASE  1 


ABC  K 


M 1WI  IHHHHHHNt  Ot  IHHHHHHHHH 


H 


F 


1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

21 

0.955567 

-0.394776 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

160 

180 

0 

32 

0.970222 

-0.261674 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

22 

0.982069 

-0.157164 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

24 

0.990826 

-0.080049 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

25 

0.996641 

-0.029226 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

26 

0.999579 

-0.003656 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

27 

0.999728 

-0.002363 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

28 

0.997192 

-0.024419 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

29 

0.992094 

-0.068947 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

40 

0.984566 

-0.135104 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

41 

0.974756 

-0.222079 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

42 

0.962821 

-0.329093 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

41 

0.948922 

-0.455388 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

44 

0.922220 

-0.600225 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

45 

0.915917 

-0.762881 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

46 

0.897156 

-0.942645 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

47 

0.877121 

-1.138807 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

48 

0.855985 

-1.350672 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

49 

0.823918 

-1.577533 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

50 

0.811082 

-1.818697 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

81 

0.261312 

-11.656804 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

82 

0.255355 

-11.857100 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

82 

0.250133 

-12.036570 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

84 

0.245633 

-12.194245 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

85 

0.241845 

-12.329250 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

86 

0.238759 

-12.403822 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

87 

0.236367 

-12.528270 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

88 

0.234663 

-12.591125 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

89 

0.233642 

-12.628973 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

90 

0.233302 

-12.641622 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

91 

0.233642 

-12.628973 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

92 

0.234663 

-12.591125 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

92 

0.236367 

-12.528270 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

94 

0.238759 

-12.440822 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

95 

0.241845 

-12.329267 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

96 

0.245633 

-12.194263 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

97 

0.250133 

-12.036578 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

«8 

0.255355 

-11.857116 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

99 

0.261311 

-11.656837 

1.0 

1.0 

1.0 

10.55 

0.50 

0.50 

0.50 

0.50 

0 

180 

180 

0 

100 

0.268014 

-11.436835 
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